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ABSTRACT 



Context. Giant cluster radio relics are thought to form at shock fronts in the course of collisions between galaxy clusters. Via processes 
that are still poorly understood, these shocks accelerate or re-accelerate cosmic-ray electrons and might amplify magnetic fields. The 
best object to study this phenomenon is the galaxy cluster CIZA J2242. 8+5301 as it shows the most undisturbed relic. By means of 
Giant Metrewave Radio Telescope (GMRT) and Westerbork Synthesis Radio Telescope (WSRT) data at seven frequencies spanning 
from 153 MHz to 2272 MHz, we study the synchrotron emission in this cluster. 

Aims. We aim at distinguishing between theoretical injection and acceleration models proposed for the formation of radio relics. We 
also study the head-tail radio sources to reveal the interplay between the merger and the cluster galaxies. 

Methods. We produced spectral index, curvature maps and radio colour-colour plots and compared our data with predictions from 
models. 

Results. We present one of the deepest 153 MHz maps of a cluster ever produced, reaching a noise level of 1.5 mjy beam -1 . We 
derive integrated spectra for four relics in the cluster, discovering extremely steep spectrum diffuse emission concentrated in multiple 
patches. We find a possible radio phoenix embedded in the relic to the south of the cluster. The spectral index of the northern relic 
retains signs of steepening from the front towards the back of the shock also at the radio frequencies below 600 MHz. The spectral 
curvature in the same relic also increases in the downstream area. The data is consistent with the Komissarov-Gubanov injection 
models, meaning that the emission we observe is produced by a single burst of spectrally-aged accelerated radio electrons. 

Key words. Acceleration of particles - Radio Continuum: galaxies - Galaxies: Clusters: individual : CIZA J2242. 8+5301- Galaxies: 
clusters: intracluster medium - Cosmology: large-scale structure of Universe 



1. Introduction 

Galaxy clusters, the most massive gravitationally-bound objects 
in the Universe, have most of their baryonic mass in the form of 
hot intra-cluster gas visible in the X-rays. In the context of hi- 
erarchical structure formation, clusters mainly grow by mergers 
with other clusters and galaxy groups, events which release into 
the ICM enormous amounts of energy up to the order of 10 64 
erg on time scales of 1-2 Gyr (e.g., Sarazin 2002). In these ex- 
treme environments, radio emission can be found in the form of 
disturbed radio galaxies, relics and haloes. 

The intra-cluster medium (ICM) interacts with the radio 
galaxies travelling at high speeds by shaping their radio jets into 
a head-tail morphology (Mile y et al.|1972| l. A spectral index gra- 
dient across the tail is expected as the electrons in the ram pres- 



et al. ( 200 1 1> and |Petrosian (2001 1 have proposed a turbulent 



re-acceleration mechanism, where cluster shocks induce tur- 
bulence, which in turn re-accelerates fossil relativistic parti- 
cles. Other studies propose that the emission comes from sec- 
ondary electrons injected by proton-proton collisions (e.g., Blasi 
& Colafrancesco|1999||Dolag & Ensslin|2000) . 



sure stripped lobes lose energy via synchrotron emission ( jO'Dea 
|et al.|1987> . 

Haloes are centrally-located Mpc-wide diffuse, unpolarised, 
steep-spectrum objects, that follow the spatial distribution of 
the ICM as seen in the X-ray. To explain their origin, Brunetti 



The main focus of this paper are relics, extended, Mpc-wide, 
diffuse, polarised radio emission in the form of arc-like, fila- 
mentary structures at the periphery of merging clusters. Their 
study can unravel not only how mergers affect structure forma- 
tion and the evolution of clusters, but also detailed physics of the 
ICM and particle acceleration mechanisms in a cosmic frame- 
work. The origin of Mpc-scale relic emission has been debated, 
but observations suggest they are exclusively found in merging 
clusters with disturbed X-ray morphology, displaced X-ray peak 
from the galaxy distribution (e.g., |Ferrariet a l. 2008 ; F eretti et al~| 
1998] > suggested that the merger produces 
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cluster wide shock waves that travel through the ICM and accel- 
erate thermal particles through the diffusive shock acceleration 
mechanism (DSA; e.g., |Drury|1983| l. Another mechanism that 
has been proposed is shock re-acceleration of relativistic fos- 
sil electrons (Mark evTtch et al.|2005| l, which is similar to DSA, 
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but instead of accelerating particles from the thermal pool, it as- 
sumes pre-accelerated mildly relativistic electrons that originate 
from, e.g., past radio galaxy activity. Both these models predict 
a connection between the injection and integrated spectral index 
and increasing spectral index and curvature into the post shock 
area. From an observational point of view, acceleration and re- 
acceleration are probably indistinguishable since the resulting 



spectra are similar (Ma rkevitch et al.||2005 Giacintucci et al. 
2008 ). From an observational point of view, relics are perme 



ated by jt/G-level magnetic fields and emit synchrotron radiation 
with a spectral index between -0.8 and -1.6, (F v cc y ff ) and 
curved radio spectra due to synchrotron and inverse Compton 
losses (e.g., |Ferrari et al. 2008; Fe retti et al.|2012| l. Models sug- 
gest that the turbulence injected by the travelling shock in the 
downstream area produces these magnetic fields, which are then 
amplified through shock compression. Simulations by Iapichino 



& Briiggen] ( |2012[ ) predict strong magnetic fields aligned with 



the shock surface of the order of 6 fiG at 0.5/? vir , which are 



in agreement with observations (e.g., Bonafede et al. 2009 
Finoguenov et al. 2010 van Weeren et al. 201 la[). The elec 



trons are accelerated to an initial power-law energy distribution 
spectrum. They emit at low radio frequencies and their emis- 
sion spectrum is directly connected to the shock parameters, 
such as the Mach number. Synchrotron and inverse Compton 
losses affect the high energy electrons more than the low en- 
ergy ones. For an initial power-law distribution of electron en- 
ergies spectral ageing causes a cutoff at the high frequency part 
of the emission spectrum. At lower frequencies the spectrum re- 



mains a power-law (Rybicki & Lightman 1979). The shape of 
the high-frequency fall-off is determined by physical processes 
such as the energy injection to the electrons and the pitch an- 
gle to the magnetic field. Both the Kardashev-Pacholczik (KP; 



Kardashev 1962 



Pacholczyk 1970 1 and the Jaffe-Perola (JP; 



Jaffe & Perola|1973[ l modelsassume a single-shot particle injec 
tion. The former considers the pitch angle to be constant in time. 
This leads to a power law fall-off with index 4/3a , i n j - 1, where 
ai n j is the injection index. By contrast, the latter assumes a con- 
tinuous isotropisation of the angle on timescales much shorter 
than the the time after which the radio emission diminishes sig- 
nificantly (10 - 100 Myr). This leads to a faster, exponential 
fall-off. More complicated models include the continuous in- 
jection model (CI; Pacholczyk| 1 970} , where fresh electrons are 
steadily injected and the Komissarov-Gubanov model (KGJP; 
Komissarov & Gubanov||1994| l which extends the JP model to 
include a finite period of freshly supplied electrons. 

To date, only ~50 examples of relic emission are known 
(e.g., |Nuza et aL[|2012| |Feretti et aL[|2012|). Roughly ten 



of these systems contain two relics symmetrically positioned 
with respect to the cluster centre, e.g. A 1240 and A2345 
d Bonafede et~alj|2009]l A3376 ( Bag chi et aL] [2006]l, A3667 
( Rottgering et al.||1997|>, PLCKG287.0 (|Bagchi et al.| |20TT), 



ZwCl 2341.1+0000 (|van Weeren et al . 2009), Zw Cl 0008+52 
(|van Weeren et al.||2011a|), MACS J1752.0+4440 ( |van Weeren 
et al.|2012a[|Bonafede et al.|2012|i, RXCJ 1 3 14.4-25 15 ( |Ventur i 
et al.|2007[ ), 0217+70 ( |Brown et al.|20TTj l. These are likely pro 



duced as a result of a major merger between two clusters of sim- 
ilar mass (e.g., Roettig er et al.|1999"t|van Weeren et al.|201 lb[ ). 

CIZA J2242.8+5301 (from here on, C2242) is an extraordi- 
nary example of a merging cluster hosting a Mpc-wide double 
relics ( |van Weeren et al.|2010] >. At a redshift of z = 0.1921, it is 
a luminous X-ray cluster with Lx = 6.8 x 10 44 erg s , measured 
between 0.1 and 2.4 keV ( |Kocevski et al.|20 07 ), and is marked 
by signatures of a major merger event. Analysis of recent XMM- 
Newton imaging by |Ogrean et al.| ( |2013| l confirm the merging 



nature of the cluster with X-ray morphology elongated along 
the north-south direction, significant substructure and hints for 
a shock front at the location of the relics. Because of the limited 
sensitivity in the X-rays, they are unable to to detect a joint den- 
sity and temperature jump. lAkamatsu & Kawahara (2013 1, via 



Suzaku data, were able to detect an X-ray density jump at the 
location of the northern relic, consistent with a Mach number 
M = 3.15 ± 0.52+°?°. The two relics are located at 1.3 Mpc 
from the cluster centre. The northern relic (RN) is 1.7 Mpc 
across and very narrow (55 kpc). The spectral index (between 
2272 MHz and 608 MHz) varies between -0.6 to -2.0 across 
the relic. This observed gradient towards the centre of the clus- 
ter is believed to be due to the spectral ageing of electrons. A 
Mach number M = 4.6^'q can be de rived from the spectral in- 
dex information ( [Landau & Lifshitz|[T959l |Sarazin||2T)0"2"] l. The 
relic has been observed to be highly polarised up to 60%. Using 
the relationship between the width of the synchrotron emitting 
region, the characteristic time scale due to spectral ageing and 
the magnetic field, B was derived to lie between 5 and 7 fiG (van 
[Weeren et al.|2010] l. The southern relic (RS) is fainter and dis- 
turbed and is connected to the northern relic by a radio halo. | van | 



Weeren et al.| (201 lb > carried out hydrodynamical simulations 



of a binary cluster merger and varied the mass ratios, impact 
parameter, viewing angle etc., as to best fit the observed param- 
eters of C2242. They concluded that the relics are seen close to 
edge-on (i < 10°) and that C2242 underwent a binary merger 
~ 1 Gyr ago, with a mass ratio of the two components between 
1.5:1 and 2.5:1. The above mentioned properties of the relic 
provide strong evidence for shock-acceleration in galaxy cluster 
shocks ( |van Weeren et al.|2010[ ). [Kang et aL] ( |2012[ ) performed 
time-dependent DSA simulations and concluded that the obser- 
vational properties of C2242 can be accounted for with both a 
shock with Mach number M = 4.5 accelerating thermal elec- 
trons and a weaker shock of M = 2.0 with a relativistic particle 
pool. 

In this paper we aim at distinguishing between theoretical in- 
jection models such as KP, JP and CI. Another process we inves- 
tigate is the (re-)acceleration of electrons emitting at long radio 
wavelengths with the goal to discriminate between shock accel- 
eration and other proposed relic formation mechanisms such as 
adiabatic compression (Ensslin & Gopal-Krishna 2001). We also 
study the cluster head-tail radio sources to reveal the interplay 
between the merger and the cluster galaxies. We therefore carry 
out a detailed radio analysis of C2242 using the Westerbork 
Synthesis Radio Telescope (WSRT) and the Giant Metrewave 
Radio Telescope (GMRT). C2242's relatively high radio surface 
brightness and large size, make it the ideal cluster for this analy- 
sis as the quality of the spectral index and curvature fits depends 
directly on the dynamic range, resolution and sensitivity of the 
radio maps employed. We present deep, high-resolution radio 
maps and the first integrated spectra, spectral index and curva- 
ture maps of four relics in C2242 to cover a frequency range of 
almost two orders of magnitude. 

The layout of the paper is as follows: in Sect. [2] we give 
an overview of the radio observations and the data reduction, 
in Sect. [3] we present radio maps and a spectral index and cur- 
vature analysis, in Sect. [4] we discuss the implications for dif- 
ferent injection and re-acceleration mechanisms. Concluding re- 
marks can be found in Sect. [5] A flat, ACDM cosmology with 
Ho = 70.5 km s _I Mpc -1 , matter density Q.m = 0.27 and dark 
energy density Qa = 0.73 is assumed (Dunkley et al. 2009}. At 
the redshift of the cluster 1 arcmin corresponds to 0.191 Mpc. 
All images are in the J2000 coordinate system. 
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Table 1. WSRT observations 



Wavelength 24.5 cm 21.5 cm 17.5 cm 13.2 cm 

Frequency range 1170-1289 MHz 1303-1460 MHz 1642-1780 MHz 2200-2344 MHz 

Observation dates Nov 17, 2009 Aug 30, Sept 4, 9, 14, 2009 Oct. 24, 2009 Sept 9, 2009 

Integration time 60 sec 60, 30, 60, 60 sec 60 sec 60 sec 

Total on-source time llh37min 36h 03min llh44min llh58min 

Usable channels per IF 56 56 56 56 

Usable bandwidth per IF 17.5 MHz 17.5 MHz 17.5 MHz 17.5 MHz 

Channel width 312.5 kHz 312.5 kHz 312.5 kHz 312.5 kHz 

Number of IFs 8 8 8 8 

Polarisation XX, YY, XY, YX XX, YY, XY, YX XX, YY, XY, YX XX, YY, XY, YX 



Table 2. WSRT image parameters 



Wavelength 


24.5 cm 


21.5 cm 


17.5 cm 


13.2 cm 


Beam size 


22.38" x 17.14" 


20.95" x 15.80" 


15.98" x 13.10" 


12.38" x 10.04" 


RMS noise (ctrms) 


51 jiJy beam -1 


43 fiJy beam -1 


34 yuJy beam -1 


41 yujy beam -1 


Dynamic range 


-4600 


-3850 


-4450 


-3250 


Mode 


multi-frequency synthesis 


mfs 


mfs 


mfs 


Weighting 


Briggs 


Briggs 


Briggs 


Briggs 


Robust 


0.5 


0.5 


0.5 


0.5 



2. Observations & Data Reduction 

In this section, we present the calibration and imaging steps per- 
formed on the WSRT and GMRT datasets. 



2.1. WSRT observations 

C2242 was observed with WSRT at four different frequencies 
(1221 MHz, 1382 MHz, 1714 MHz and 2272 MHz) in August - 
November 2009. Eight frequency bands of 20.0 MHz bandwidth 
and 64 channels were centred around each of the frequencies. 
The observational details can be found in Table Q] 

We used the Astronomical Image Processing System 
(AIPSfTand the Common Astronomy Software Applications 
(CAS Ar] package to reduce the data. Data for each of the four 
frequencies were independently flagged and calibrated. We set 
the fluxes of the calibrators 3C 286 and 3C 147 according to the 
( 1999| > extension to the Baars et al. ( 1977| l scale. 



where constant c — 68, r is the distance from the pointing center 
in degrees and v is the observing frequency in GHz. We adopt am 



Perley & Taylor 



Initial phase solutions were obtained for the centre of each IF 
and were subsequently used to compute delay phase corrections. 
Twelve edge channels were excluded due to bandpass roll-off. 
The bandpass and phase solutions derived for the usable range 
of channels were transferred to C2242. 

The science target data were split and all the IFs belonging to 
the same frequency setup (see Table TJ were combined for imag- 
ing (for imaging parameters, see Table |2j. We reached dynamic 
ranges of above 3000 measured as the ratio between the peak in 
the map and the RMS noise. We then performed three rounds of 
phase only self-calibration and one round of phase and ampli- 
tude self-calibration. 

The final image was produced using "Briggs" weighting (ro- 
bust set to 0.5, Briggs 1995 ). The images were corrected for pri- 
mary beam attenuatioij^] 



A(r) — cos (cvr) 



(1) 



1 http://aips.nrao.edu/ 

2 http://casa.nrao.edu/ 
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uncertainty of 10% in the overall flux scale (e.g., Schoenmakers 



According to the WSRT Guide to observations 



et al.| 1998 ; Rengelink et al. 1997 1, which includes errors such as 
imperfect calibration, positional errors and flux calibration errors 
for different observations. 



2.2. GMRT observations 
2.2.1. 153 MHz Observations 

The target C2242 was observed with the GMRT at 153 MHz on 
Nov 2, 2010 during a single 8.5 hour daytime observing session, 
recording visibilities over an 8.5 MHz bandwidth. C2242 was 
observed in scans of -1 hour, interleaved with 3 minute scans 
on phase calibrator J2 148+611. The observation schedule also 
included initial and final scans on flux calibrators 3C286 and 
3C 48, respectively. 

Data reduction was performed also making use of the Source 
Peeling and Atmospheric Modelling (SPAM) package (Interna 
et al. 2009), a Python-based extension to AIPS (Grei sen |2003 1 
including direction-dependent (mainly ionospheric) calibration 
and imaging. 3C48 was used as the single flux and bandpass 
calibrator, adopting a total flux of 63.4 Jy. The bandpass filter 
edges and channels with strong RFI were flagged, yielding an 
effective bandwidth of 5.5 MHz. 3C 48 was also used to deter- 
mine the phase offsets between polarisations, and for estimating 
the instrumental, antenna-based phase offsets (needed for iono- 
spheric calibration; see Interna et al. 2009). 

All calibrations derived from 3C48 were transferred to the 
target data. After clipping excessive visibility amplitudes and 
subtracting low-level RFI ( |Athreya|2009[ ), visibilities were con- 
verted to Stokes I and averaged per 6 channels for more effi- 
cient processing, while keeping enough spectral resolution to 
suppress bandwidth smearing. The target was initially phase- 
calibrated against a point source model derived from the NVSS, 
VLSS and WENSS surveys ( |Condon et aL|[T998t [Cohen et al. 
2007 Rengelink et al.|1997 1. This gives a starting model for the 
calibration and ensures fast convergence. This was followed by 
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Table 3. GMRT observations 



Frequency range 


150-156 MHz 


309-324 MHz 


323-328 MHz 


593-623 MHz 


Observation dates 


Nov 2, 2010 


Jul 2, 2010 


Jul 3, 2010 


Nov 20, 2009 


Integration time 


4 sec 


8 sec 


8 sec 


8 sec 


Total on-source time 


6h 15min 


4h 35min 


4h 35min 


8h llmin 


Total bandwidth 


8.5 MHz 


16.6 MHz 


16.6 MHz 


32 MHz 


Channels 


128 


256 


256 


512 


Usable bandwidth 


5.5 MHz 


14.7 MHz 


14.7 MHz 


30 MHz 


Usable channels 


88 


226 


226 


480 


Channel width 


62.5 kHz 


62.5 kHz 


62.5 kHz 


62.5 kHz 


Polarisation 


RR,LL 


RR,LL,RL,LR 


RR,LL,RL,LR 


RR,LL,RL,LR 



Table 4. GMRT image parameters 



Frequency 153 MHz 



323 MHz 



608 MHz 



Beam size 

c RMS 

Dyn. range 
Mode 
Weighting 
Robust 
Grid mode 
w projection 



28.4" x 23.6" 

1.5 mjy beam -1 

-650 

mfs 

Briggs 

0.5 

widefield 
150 



12.3" x 11.3" 

0.2 mJy beam" 

-1180 

mfs 

Briggs 

0.5 

widefield 
150 



6.9" x 5.1" 

24 pjy beam -1 

-2750 

mfs 

Briggs 

0.5 

widefield 
150 



rounds of (phase-only) self-calibration and wide-field imaging. 
Bright outlier sources were included in the imaging. Additional 
flagging of bad data was performed in between rounds. 

After a last round of (amplitude & phase) self-calibration 
and imaging, we performed two additional rounds of SPAM 
ionospheric calibration and imaging. This consists of "peeling" 
bright sources within the field-of-view, fitting a time and spa- 
tially varying ionosphere model to the resulting gain phase so- 
lutions, and applying this model during imaging. A peeling cy- 
cle consists of obtaining directional dependent gain solutions for 
the "peeled" source, which is now properly modelled and can 
be completely removed from the uv data (e.g., Noordam 2004). 
We peeled 16 sources in the first round and 17 sources in the 
second round of SPAM. This procedure removed a substantial 
fraction of the remaining residual artefacts around point sources 
after self-calibration. We imaged the data using the polyhedron 
method (Perley 1 989] |Corn well & PerleyfT992] l to minimise ef- 
fects of non-coplanar baselines. The final image (see Table 
made using robust 0.5 weighting (Briggs 1995), has a central 
background RMS noise level of 1.5 mJvbeam" 1 . The image was 
corrected for primary beam attenuation*] 



A , , , -3.397 , 47.192 4 -30.931 6 7.803 s , ns 
A(x) = 1 + — x 2 + x 4 + — x 6 + — pr-X 8 , (2) 



10 3 



10 7 



10 1 



10 1 



with x the distance from the pointing centre in arcmin times the 
observing frequency in GHz. 

We adopt an uncertainty in the calibration of the absolute 
flux-scale for GMRT of 10% ( |Chandra et al.|2004| >. 

2.2.2. 316 MHz and 330 MHz Observations 

C2242 was observed with the GMRT for a total of 15 hours 
on two subsequent nights, July 2 and July 3, 2010. Data were 
recorded in two slightly overlapping spectral windows centred 



at frequencies 316 MHz and 330 MHz. The 16.6 MHz band- 
width per spectral window was sampled in 256 channels in full 
polarisation mode. The observations are summarised in Table [3] 
3C 48 and 3C 147 were used as flux calibrators and 2355+498 as 
phase calibrator. 

We used CAS A to reduce the data. The two spectral windows 
were independently flagged and calibrated. The first integration 
of each scan was flagged to account for system instabilities when 
moving between set-ups. On both nights there were malfunction- 
ing antennas, resulting in 26 working antennas during the entire 
two nights of observations. The RFI in the calibrator data was 
manually removed. We set the fluxes of the primary calibrators 
according to the ( |Perley & Taylor|1999| l scale. Initial phase so- 
lutions were obtained for a small number of channels free of 
RFI close to the centre of each bandpass, which were used to 
compute delay phase corrections. Thirty edge channels were not 
considered in the calibration due to bandpass roll-off, resulting 
in 14.7 MHz of effective bandwidth per spectral window. The 
bandpass and phase solutions derived for the full usable range of 
channels were transferred to C2242. 

The science target data were split and averaged down to 
56 channels per spectral window, which were afterwards com- 
bined for imaging (for imaging parameters, see Table The 
egregious RFI was visually identified and excised. We then per- 
formed three rounds of phase only self-calibration and one round 
of phase and amplitude self-calibration. 

A CLEANed image was produced from the flagged dataset, 
whose dynamic range was limited by the presence of very bright 
sources increasing the RMS in the field. The brightest sources 
were successively removed through the "peeling" method. The 
final image was produced using "Briggs" weighting (robust set 
to 0.5, Briggs 1995[ ). Reasons for not reaching the thermal 
RMS include: residual RFI and calibration errors (Bhatnagar 



According to the GMRT User's manual 



|et al.|2008 1. The residual patterns around the positions of bright 
sources are caused by incomplete peeling due to imperfect mod- 
els of the sources, pointing errors and non-circularity of the an- 
tenna beam pattern. 

2.2.3. 608 MHz observations 

The cluster was observed at 608 MHz for over 8 hours on Nov 
20, 2009 in full polarisation (see Table HJ. 3C286 and 3C48 
served as flux calibrators, while phase calibrator 2148+611 was 
observed every hour. 

CASA and AIPS were used to reduce the data in a similar 
fashion as presented in Sect. |2.2.1 After peeling the brightest 
sources, we performed an additional step in simply subtracting 
from the uv data the contaminating sources towards the edges of 
the field. We then imaged the central part of the FOV with the 
parameters given in Table |4] 
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Fig. 1. GMRT 153 MHz image with contours drawn at 
[3, 8, 16, 32] xctrmj. The beam size at 28.41" x 23.60" is shown 
in the bottom left comer of the image. 



Fig. 3. GMRT 608 MHz image with contour lines drawn at 
[4, 8, 16, 32] x cr RMS . The beam size is 6.87" x 5.12" and shown 
at the bottom left corner of the image. 




42' ' 3(r 



Right Ascension 

Fig. 2. GMRT 323 MHz image with overlaid contours at 
[4, 8, 16, 32] xcr RM s. The beam size at 12.26" x 1 1.27" is shown 
in the bottom left comer of the image. Added source labelling. 



3. Results 

In this section we present radio continuum maps, spectral index, 
spectral curvature and colour-colour plots. 



3.1. Radio continuum maps 



van Weeren et al. (2010) presented a 1400 MHz image in their 
paper focussing on a discussion of two counter relics, pointing 



out the northern relic's narrow and elongated structure. Here, we 
discuss in detail the morphology of all the compact and diffuse 
radio sources in the field. 



3.1.1. Relics 

All of the radio maps present the two counter relics and two 
smaller areas of diffuse emission, plus a variety of radio sources 
(for labelling see Figs. [T]and|4]i. The sizes reported in the fol- 
lowing text are measurements in the high-resolution 608 MHz 
image within 5cr RMS levels and represent true sizes convolved 
with the synthesised beam at this particular frequency. 

The northern relic (labelled RN), maintains its arc-like shape 
and impressive size (1.8 Mpc x 150 kpc) throughout more than 
one order of magnitude in frequency. The thickness of the relic 
in the north-south direction increases with decreasing frequency 
(~160 kpc), as expected for a source with a strong spectral index 
gradient across it. The tightness of the contours reveals the north- 
ern boundary to be extremely well defined and sharp (see Fig. 3}. 
Towards the inner edge, the surface brightness drops smoothly 
and slowly. Towards the east of RN, source H has a much higher 
surface brightness than the relic at all frequencies and presents a 
typical tailed radio source morphology towards the north. 

The southern relic (RS) does not posses the same well de- 
fined structure as RN, but generally follows a bow geometry. 
Its projected size is 590 kpc x 310 kpc. In the 2272 MHz map 
(Fig. [7]) the relic is very faint, because of its steep spectrum. The 
WSRT and 608 MHz maps reveal a tail of fainter, diffuse emis- 
sion extending towards the south that is not visible in the 153 
and 323 MHz images (Figs.[T]and|2|. At first glance this is puz- 
zling if we do not take into consideration the noise levels of the 
GMRT measurements, which are an order of magnitude higher 
than the other maps. The spectrum of this extension should be 
steeper than a - - 1 .4 for it to be detected in the low frequency 
maps. 
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Fig. 4. Left: WSRT 1221 MHz image with contour lines drawn at 
[4, 8, 16, 32] x ctrms. The beam size at 22.38" x 17.14" is shown 
in the bottom left comer of the image. Added source labelling. 
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Fig. 5. WSRT 1382 MHz image with contour levels at 
[4, 8, 16, 32] x ctrms. The beam size at 20.95" x 15.80" is shown 
in the bottom left comer of the image. 



To the north of RS, the 1221 MHz map (Fig. |4]i reveals a 
patch of diffuse emission labelled as source K. The patch gets 
connected to RS in the 1221 and 1382 MHz maps (Fig. [| and [5}. 

Diffuse source J is located at the west of RS and covers an 
area of 260 kpc x 350 kpc. In the lower frequency maps (below 
1382 MHz), source J and RS become connected into a single 
area of diffuse emission with flux concentrated in two patches. 
The distinction between source J and RS will become evident in 
the spectral index maps (see Sect. |3 .2.2) . 




Right Ascension 



Fig. 6. WSRT 1714 MHz image with contours at [4, 8, 16, 32] x 
0"rms- The beam size at 15.98" x 13.10" is shown in the bottom 
left corner of the image. 
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Fig. 7. WSRT 2272 MHz image with contour lines drawn at 
[4, 8, 16, 32] X 0-rms- The beam size at 12.38" x 10.04" is shown 
in the bottom left corner of the image. 



Towards the east of RN, the maps reveal another smaller 
(350 kpc x 500 kpc) arc-like area of extended, diffuse emission 
(Rl). 

In all of the maps, a patch of extended emission (R2) 
is detected immediately south of the northern relic, with a 
size of 380 kpc x 290 kpc. In the high frequency maps 
(1221, 1382 and 1714 MHz) the emission has two peaks, but the 
323 MHz map reveals that the western peak is actually a separate 
point source. This source, embedded in the relic emission, has a 
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Table 5. Integrated radio spectra of the relics and steep spectrum source G 



Source 


RN 


RS Rl 


R2 J G 


a 


-1.06 ±0.04 


-1.29 ±0.04 -0.74 ±0.07 


-0.90 ± 0.06 -1.53 ± 0.04 -1.77 ± 0.05 




Frequency [MHz] 



Fig. 8. Integrated radio spectrum of RN. 




Frequency [MHz] 



Fig. 9. Integrated radio spectrum of RS. 



0.1 00 




Frequency [MHz] 



Fig. 10. Integrated radio spectrum of Rl. 




Frequency [MHz] 

Fig. 11. Integrated radio spectrum of R2. 

1 .000 E 




Frequency [MHz] 

Fig. 12. Integrated radio spectrum of source J. 

0.1000E 




Frequency [MHz] 



Fig. 13. Integrated radio spectrum of source G. 
3.1.2. Radio galaxies 



steeper, brighter spectrum than the surrounding diffuse emission Within the sensitivity limits of our observations we detect five ra- 

which, at some frequencies, is below the noise levels. Towards dio sources in the field with head-tail morphology. The 608 MHz 

the west, immediately south of RN, source I is an elongated image (Fig. [3j reveals object E to be a classical twin-tailed radio 

patch of emission, which is visible only in the 608 MHz, the source with highly bent lobes. Interesting is also source G which 

1221 MHz and 1382 MHz images (Figs.[3]|4]and|5]l. is at the noise level in the 2272 MHz map, but a strong detection 
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Fig. 14. Three-frequency spectral index map between 153 MHz and 608 MHz. Contours from 1221 MHz are overlaid at 
[4,8, 16,32] xo-rms level. 



at the low frequencies. This means the source has an extremely 
steep spectrum. It has to be noted that source A has been peeled 
in the 153 MHz and 608 MHz. Since the source could not be per- 
fectly removed from the uv data, the source itself is not visible, 
but there is some residual side lobe pattern in the images. 

3.2. Spectral index and spectral curvature 

For the purposes of spectral index and curvature mapping, it is 
important to have consistent flux scales for all the seven fre- 
quencies. Despite the fact that the same flux standard was used 
throughout the calibration process, issues with the transfer of the 
flux scale from the calibrator to the target and differences in uv 
coverage result in the recovery of flux on scales that vary be- 
tween observations and telescopes. 

In order to account for this, we produced new radio maps 
with a common, minimum uv cut at Q.2kA and uniform weight- 
ing to ensure recovery of flux on similar spatial scales. This uv 
cut accounts for the missing short uv spacings and causes the 
largest detectable scale to be ~ 17'. The resolution of these maps 



is, in ascending order of frequency: 16.2" x 13.8", pa = 50.0°; 
10.0" x 8.0", pa = 0°; 4.2" x 3.4", pa = 60.0°; 16.9" x 12.9", 
pa = -1°; 14.9" x 11.0", pa = -1°; 12.0" x 9.1", pa = 0.0°; 
9.5" x 7.1", pa - 0.0°. We then convolved all the maps to the 
lowest resolution and regridded/aligned them with respect to the 
1221 MHz image. The final effective resolution is: 18.0"xl4.8", 
pa = 20.0°. 

The uncertainty in the flux scale for both the WSRT and 
GMRT observations was considered 10%. We base this value 

1998] l 



on the studies of Schoenmakers et al. 



Rengelink et al. 



([1997} and |Chandra et al.| ( |2004) >. To~ensure that the uncertainty 
is not underestimated, we studied the spectrum of seven compact 
sources in the field, using our measurements together with values 
from the NVS S (|Condon et al.|1998| , WENSS ( |Rengelink etaL 



1997b, VLSSI Cohen eta 



1996), 4C (Gower et al. 



2007| >, Texas 365 MHz (|Douglas et al 
1967l land 7C ( |Vessey & Green||1998) 



catalogues, where available. We derived linear fits to the mea- 
surements as function of frequency and measured the dispersion 
of the points with respect to the best fit. We concluded that our 
measurements were all within 10% of the fitted lines. 
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Fig. 15. Seven-frequency spectral index map between 153 MHz and 2272 MHz. Contours from 1221 MHz are overlaid at 
[4,8,16, 32JXO-RMS level. 



3.2.1. Integrated spectra 

We produced integrated radio spectra for the four relics RN, RS, 
Rl and R2 using these maps with the common uv cut and same 
resolution. We measured the fluxes in fixed-size apertures across 
the seven frequencies. The uncertainty in the flux was computed 
as the RMS noise multiplied with the number of beams A^eams 
contained in the respective area and the 10% flux uncertainty, 
added in quadrature: 



AF= V(^RMs) 2 JVbea mS + 0.01F 2 . 



(3) 



Figure [8] shows the flux measurements for RN, with a linear 
fit. Although the data is well described by a single power law, 
the spectrum falls off at 2272 MHz and becomes more curved. 
Figures|9 10 and [TT|present the linear fits for RS, Rl and R2, re- 
spectively. Rl is the only relic not well described by the straight 
spectrum. Diffuse source J is fitted with a power law with slope 
-1.55. The details of the fits for the relics are given in Table [5] 

We also produced a linear fit to the spectrum of source G 
using six frequencies. Since it is not detected in the 2272 MHz 



image, we plotted, for reference, a 3cr RM g upper limit value at 
this frequency (see Fig.[T3). 



3.2.2. Spectral index maps 

For the spectral index maps, we blanked all pixels with values 
lower than 2cr RM s and fitted linear functions to the data pixel 
by pixel. Errors were computed by adding in quadrature the flux 
uncertainty error of 10% and the cr RMS in each map. 



Figure 14 shows the low frequency GMRT spectral index 



map between 153 and 608 MHz. A spectral index gradient across 
the northern relic is immediately evident, varying from an av- 
erage value of -0.6 at the outer to -2.5 at the inner edge, as 
was first observed at high radio frequencies by |van Weeren et aT] 
(2010). The gradient is smooth and consistently perpendicular to 
the relic from its east to its west tip. The beam is sampled only a 
few times over relic's thickness and, at the northern side, there is 
an abrupt drop in surface brightness where edge effects become 
important. This causes a sharp decrease in the spectral index in 
parts of the northern side of the relic. There is some spectral 
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Fig. 16. Seven frequency spectral curvature map between 153 MHz and 2272 MHz. Contours from 1221 MHz are overlaid at 
[4,8, 16,32] xo-rms level. 



index steepening across relic Rl, while R2 has a steep spectral 
index {a < -1.3), but no strong gradient. A spectral index gradi- 
ent is also noticeable across the head-tail sources B, C, D, E and 

F. The nuclei of these sources have a 0.5, whereas the tails 

steepen to values below -2.0. The most dramatic steepening can 
be observed across source F, whose tail reaches spectral index 
values of -3.5. Source G, which was below the noise level in 

the 2272 MHz map, has an ultra steep spectrum with a 1 

The nuclei of tailed radio sources have a typical spectral indices 

f!994l l 



of ~ -0.7 with tails steepening to -2. Rottgering et al. 



classify all sources with integrated spectral index a < -1.0 as 
ultra steep spectrum sources. The point source embedded in the 
western part of R2 is now properly resolved. The spectral in- 
dex is much flatter than the relic's and fairly constant across the 
source. 



In Figure 15 



we present seven frequency spectral index 
maps with a frequency coverage from 153 MHz to 2272 MHz. 
As before, the gradient across the northern relic is clearly visible. 
Although in the low-frequency radio maps, the southern relic 
and source J appear as a single object, in the spectral index map 



they are split in two areas with different properties. RS, while not 
possessing the same orderly morphology as RN, still displays a 
gradient with steepening spectrum (-3.0 < a < -1.5) towards 
the cluster centre. Source J has a much steep spectral index 
(-2.0 < a < -1.7), while the RS is flatter (-1.2 < a < -0.6). 



g 3.2.3. Spectral curvature map 

We produced a spectral curvature map using all of the seven 
maps available convolved to the same resolution. The wide fre- 
quency coverage enables us to better constrain the curvature and 



minimise errors. We follow the definition of Leahy & Roger 



(1998 ) for three frequencies: 



C 



-a,' + al; 



(4) 



where v\ is the lowest of the three frequencies, V2 is the center 
one and V3 is the highest. In this convention, the curvature C for 
a normal, convex spectrum (e.g. spectral ageing) is negative. As 
we were unable to find a consistent definition of curvature when 
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Frequency [MHz] 
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Fig. 17. Normalised spectrum with linear fits for the regions selected based on spectral index. The flux was multiplied by the 
frequency on the vertical axis to better differentiate between the different regions. The data points were slightly shifted along the 
horizontal axis for clarity. 



more than three frequencies are employed, we define spectral 
curvature in the following way: 



C 



1650 
380 



-£*low v + ffhigh,. 



(5) 



where the low frequency spectral index was computed using 
the three GMRT frequencies centred at 380 MHz and the high 
frequency spectral index was obtained from the fit to the four 
WSRT measurements, centred at 1650 MHz. Thus the curvature 
depends on which frequencies one uses for the spectral indices. 
The two fits also produced standard errors of the spectral index 
(Aai ow „ and Aarhi g i,_v) which were added in quadrature to obtain 
the uncertainty in the curvature estimate: 



AC 



1650 
380 



(Aaio 



„) 2 + 



(Aahigh v ) 



(6) 



The results of a pixel by pixel curvature is presented in 
Fig. [16] We blanked all pixels with SNR smaller than 3. As ex- 
pected, the relics and the tails of the radio galaxies have a convex 
spectrum, i.e. a negative curvature parameter. 

RN is marked by a gradient of increasing curvature from 
north to south. At the front of the shock we have C — 0, which re- 
inforces the results of the separate high and low frequency spec- 
tral index maps. The curvature increases into the downstream 
areas to values of -1.5. There are also small scale variations 
along the length of the source with a size of the order of a 
beam (~64 kpc). We treat these in more detail in the Discussion 
(Sect |) 

Rl and R2 do not show any extreme curvature, while the di- 
vision between source J and RS is still noticeable from their dif- 
ferent curvature properties. The spectrum of RS is considerably 



less curved than J which reaches C ~ -1.5. There is smooth 
gradient over the tail of source E, which dips to -0.7. A pro- 
nounced spectral curvature is also visible in the tail of source F: 
values range from -1.0 to -1.4. We measure source G separately 
and it is also extremely curved with C 1.25. 



3.3. Colour-colour plots 

The spectral index and curvature maps were produced on a pixel 
by pixel basis. In order to increase the signal to noise ratio 
(SNR), we perform a spectral curvature analysis for RN on re- 
gion with similar spectral properties based on the method of van 
eeren et al. ( 20 1 2b| > . For this, we bin pixels in seven spectral 
index groups ranging from -0.75 ... - 1 .35, based on their seven- 



frequency a in the spectral index map. The step size was chosen 
as 0.1, to gather enough pixels (> 800, pixel size is 1" x 1") to 
get high SNR in each bin, but to avoid too much mixing of elec- 
trons from different age populations. As compared to a pixel by 
pixel analysis, the binning increases the SNR with respect to the 
RMS at least a factor of 200. 

The northern, outer edge of the relic is the front of the shock. 
The spectral index increases with distance into the back of the 
shock, towards the inner edge. Therefore, the spectral index cri- 
terion traces regions in a shell-like pattern from the shock into 
the post-shock area. 

The spectral index selection criterion produces regions that 
vary in size, which makes it difficult to produce directly com- 
parable spectra for each of these areas of the relic. In order to 
account for the difference in surface area, after we sum up the 
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Fig. 18. Curvature as function of seven-frequency spectral index 
for different regions in the northern relic. 

flux fj of all of the pixels in each area, we normalise by the num- 
ber of pixels to get an average flux per pixel: 



N pixels 



(7) 



where F is the average flux for the region. The flux uncertainty 
derived in Sect. |3.1.T| is then normalised by the number of pixels 
in the respective area to get the standard error of the normalised 
flux F. We fit second order polynomials to the normalised fluxes 
as a function of frequency. We then use the best-fit parameters 
to predict the flux at our reference frequencies. We thus have 
spectra predicted at seven frequencies for all the eight regions. In 
this way we are improving the SNR, by using all seven available 
flux measurements jointly to predict the flux at each of the seven 
frequencies. In the next analyses, we use these predicted fluxes 
for each of the a-selected regions, rather than the measurements 
themselves. 

The spectra for the different regions of RN are plotted in 
Fig. 17 where the predicted flux was multiplied by the frequency 



on the y-axis to emphasise the differences between the regions. 
Linear fits were drawn through the points for reference. The 
colours from black to red represent areas selected based on the 
pixels increasing spectral index. We compare the flux points with 
the fit in order to evaluate how strongly curved the RN spectrum 
is. The slope of the linear fit corresponds to the spectral index 
that describes the pixels summed up within an area. The areas 
with flatter spectra are well described by a power law, while the 
spectra become increasingly curved with steepening a and devi- 
ate from the linear fit. 

The curvature was then computed using the formula defined 
in equation [5] with the same reference frequencies as before. In 
order to better visualise the dependence of the curvature on the 
spectral index, we have plotted these two quantities against each 
other in Fig. 18 The plot shows that the spectrum becomes more 



curved further from the shock. The dependence of the curvature 
on the spectral index fitted through the seven frequencies can be 
described by a linear function of the form C^ = (0.7 + 0.1) ■ 
a™ 4 + (0.3 ±0.1). 



3.3.1. Injection in the northern relic 

To test how well the RN data is described by some of the estab- 
lished injection models (see Sect.[T]i, we produced colour-colour 
plots (Katz-Stone et al. 1993 ), in which the high-frequency spec- 
tral index for multiple positions in the radio source is plot- 
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Fig. 19. Colour-colour diagram showing RN data points in black 
dots. Different injection models have been overplotted for refer- 
ence. 



ted against the low-frequency spectral index. Colour-colour di- 
agrams have the advantage of emphasising spectral curvature 
and displaying data for all areas of the source in an empiri- 
cal way. They can reveal trends which can be overlooked when 
fitting physical models directly to the data. They represent an 
easy way to visualise the models and to trace back the data 
to injection properties. This is possible because the shapes of 
spectral models are conserved under adiabatic and magnetic 
field changes and radiation losses. For this purpose, we used 
the same regions defined in Sect. |3.3| and fitted the low fre- 
quency (between 153 and 608 MHz) and high frequency (be- 
tween 1221 and 2272 MHz) spectral index using the predicted 
fluxes from the second order fit to the data. For reference, 
we overplotted JR KP, KGJP and CI models described in the 
Introduction (Sect. [TJ with injection spectral index of -0.6 and 
-0.7, in order to match the injection index derived from the spec- 
tral index maps. The intersection of the traced back data to the 



608 
r 153 



a??!? line gives an injection spectral index between -0.6 



and -0.7, which is consistent with the one derived from the spec- 
tral index map. 



4. Discussion 

The northern relic in merging cluster CIZA J2242. 8+5301 has 
been one of the best studied objects of its type owing to the size, 
regular morphology and high surface brightness. |van Weeren 
et al. ( |2010| l have performed a high radio frequency analysis fo- 
cussed on the northern relic, discovered spectral index steepen- 
ing and aligned magnetic field vecto rs and derived a Mach num- 
ber of M 



4 6+1.3 
* M -os- 



van Weeren et al. 



(2011b) 



Simulations of 

managed to reproduce the morphology and spectral index trends 
of RN within a head-on collision of two similar-mass clusters 
in the plane of the sky. The merger produces the main counter- 
relics via opposite travelling shock waves. 
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The new low-frequency radio data combined with the exist- 
ing GHz measurements provide an excellent opportunity to ex- 
tend previous research to a cluster-wide analysis. This enables 
us to answer some outstanding questions about spectral curva- 
ture and injection/acceleration mechanism. We will discuss the 
morphology of the sources across seven radio maps to provide 
clues about the nature and origin of the sources. We will interpret 
the spectral index and curvature map to fix the physical prescrip- 
tions for the four diffuse sources in the cluster, such as the Mach 
number and formation mechanism. The colour-colour analysis 
will focus on determining the injection mechanism responsible 
for accelerating the electrons within the northern relic. 

4.1. Northern relic 
4.1.1. Morphology 

The northern relic possesses the morphological and spectral 
characteristics of shock wave induced emission. It maintains its 
arc -like shape over 1 Mpc and almost two orders of magnitude in 
frequency. Source H toward the east is most probably a separate 
physical system, whose position coincides with the projected lo- 
cation of the shock. 



4.1.2. Spectral analysis 



In the acceleration scenario (En sslin et al.||1998) , at the front 
of the travelling shock we expect a straight and fiat spectrum, 
as the entire electron population is similarly accelerated. The 
spectral index a at the outer edge of the RN, where accelera- 
tion is actively happening, has a value ~ -0.6 from 153 MHz 
to 2272 MHz. The injection spectral index between -0.6 and 
-0.7 derived from the spectral index maps matches the colour- 
colour analysis. Moreover, the spectra for regions selected based 

on spectral index in Fig.fTTlcan all be traced back to 0.60 

at the low frequency encTT where radiation losses have not af- 
fected the spectral shape. The integrated spectrum of RN is well 
described by a power law with -1.06 + 0.05 slope and does 
not show any flattening or turnover at low frequencies, which 
is expected from shock acceleration theory. This is confirmed 
by the fact that the or™ derived from the low-frequency, high- 
frequency and seven-frequency map derived are consistent with 
each other at a value of 0.60. is 0.5 flatter than the inte- 
grated spectral index which is expected for a simple shock model 
(e.g., Miniati 2002, Bagchi et al. 2002). The acceleration mech- 
anism also predicts increasing spectral curvature with depth into 
the downstream area. We observe large scale trends of increas- 
ing curvature in our maps (Fig. 16 1 as well as in our colour- 
colour plot (Fig. fig). The curvature at the front of the shock is 



consistent with 0, which indicates that all electrons emitting be- 
tween 153 MHz and 2272 MHz are currently accelerated. The 
curvatures increases to = -ai OWr + a/,i g h v = -1.5 in the 

downstream area, showing a curved spectrum marked by energy 
losses. 



4.1.3. Relic substructure 

We found the significant curvature substructure along the length 
of the northern relic particularly puzzling. Are the variations due 
to real physical conditions or just noise peaks? In order to inves- 
tigate the origin of these so-called "ripples" we placed circular 
bins across the relic, with a size matching the synthesised beam 
(radius~32 kpc). Their positions were chosen such that they fol- 
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Fig. 20. Top: areas selected for measuring the spectral curvature 
and flux. Bottom: Fluxes with error bars measured within the 
circular regions for seven frequencies 153 MHz to 2272 MHz. 
The x-axis measured distance along the relic measured in kpc 
from the first region to the east, towards west. The curves in 
solid lines represent increasing frequency from top to bottom. 
Average fluxes between 323 MHz and 1714 MHz are overplot- 
ted in a dashed red line. The curvature within the same regions 
is plotted in a dashed black line. The curvature values were 
shifted along the y-axis for plotting reasons and stretched for 
clarity. The colours in the solid line at the top of the plot indi- 
cate the level of confidence for the flux trends: red correlates to a 
noise dominated measurement, while green represents trustwor- 
thy measurements. 



low the surface brightness distribution contours along the relic 
(see top map in Fig. 120). 

In the bottom plot in Fig. [20] we plotted the flux measured 
within these circular regions. It is immediately obvious that 
the variations within the inner five frequencies (323 MHz to 
1714 MHz) correlate well and show similar trends. This sug- 
gests there is real substructure associated with the shock. The 
radio maps at the bounding frequencies 153 MHz and 2272 MHz 
were heavily affected by the uv cut imposed for spectral analysis 
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purposes. This leads to a reduction in relic brightness with re- 
spect to the background noise. We suspect spurious peaks in the 
brightness distribution across there two frequencies affected the 
curvature the most. To test this, we look at the curvature varia- 
tions along the relic, fitted using the fluxes measured in the same 
circular regions. We detect suspicious maxima and minima at 
distances ~450 kpc, ~900 kpc and ~ 1 150 kpc, which can be 
traced directly to noise peaks in the 153 MHz and 2272 MHz 
maps. We therefore suggest that the low curvature regions at 
RA = 53°09'15", DE = 22 /, 43 m 02* and RA = 53° 10' 00", 
DE = 22 /! 42 m 39 I and the curvature trough at RA = 53°09'45", 
DE = 22 /l 42"'48 s are probably noise effects. We have indicated 
the areas with poor signal to noise and low degree of confidence 
in red at the top of Fig. 



20 



We notice the variations are stable on scales of roughly 
100 kpc by following the average flux trends and the curvature 
variations. These are marked in green in the "confidence bar" 
in Fig. 



20 



Focussing on these correlated flux trends, we spec- 
ulate the variations indicate complexity in the shock structure 
on similar scales. As shown by simulations (|van Weeren et al.' 



|2011b[ ), the shock has travelled for hundreds of Myr since the 
core passage between the merging clusters. Breaking of shock 
surfaces along their path could result in brightness variations 
along relics when projected onto the sky (Vazza et al. 2012| l. 
Since the radio emission is most likely produced by electrons 
accelerated by the shock within a magnetic environment, vari- 
ations in ICM density and temperature result in wiggles in the 



synchrotron emission/surface brightness along the relic (Hoeft 
et al.|20lT i. Simulations by Kang et al. (20121 also suggest that 



the surface of the shock that creates the northern relic should be 
highly elongated and that it would not be trivial to induce such a 
structure in cluster mergers. They argue that the relic might con- 
sist of a number of substructures, which is consistent with our 
observations. 



4.1.4. Acceleration of emitting particles 

The accepted model for acceleration of particles in the giant ra- 
dio relic RN is (re-)acceleration (Ensslin et al. 1998 ; Mark evitch| 
et al.|200"5| l. DSA and re-acceleration require the presence of a 
large scale shock wave and predict a single power law spectrum 
at the low frequencies, with steepening at the high frequencies 
due to energy loses into the post-shock area. There are two pos- 
sible sources for the shock: accretion shocks at the interface be- 
tween infalling material and an already formed cluster structure 
or merger shock that occur in major galaxy cluster mergers. The 
former predicts high Mach numbers (> 4) and occur at large dis- 
tance from the cluster centre in low density environments, while 



the latter predicts average Mach numbers (2 < M < 4) ( Bykov 
|et al.|20 08 and references therein). 

As mentioned in Sect. [T] after the passage of the shock, the 
accelerated electrons start losing energy via inverse Compton 
and synchrotron processes, which affect the high energy parti- 
cles first. For the physical parameters of RN, using the formal- 
ism from Hoeft & Briiggen ( 2007| l, we can derive an electron 
cooling time for the northern relic of a few 10 Myr, which peaks 
for a magnetic field of a few //G. By multiplying this time scale 
with the shock velocity, we obtain at the distance at which the 
electrons start to cool significantly: ~ 50 kpc, which means in 
the northern relic we can resolve the downstream cooling. 

The age of the electron population increases away from the 
shock, such that lower energy electrons also get affected. The 
decreasing spectral curvature at the back of the shock is to be ex- 
pected as the cut-off frequency where the electrons start losing 



energy is shifted towards lower frequencies. The spectral cur- 
vature map (see Fig. 16 1 provides strong evidence that shock 
acceleration is the responsible physical process for RN, with no 
spectral curvature at the shock front and increasing curvature in 
the post-shock area. In the context of shock acceleration, the in- 
jection spectral index is conn ected to the Mach number M of 

M 2 +3 



the shock via a 



Weeren et al. 



inj 



2M 2 -2 



(Blandford & Eichler 



1987) 



van 



4.6+JI using 



(2010) have obtained a value of M 
high frequency data only. Our results confirm their findings into 
the low frequency regime. Assuming an injection spectral index 
of cfinj = -0.60 ± 0.05, as measured from the spectral index 
maps and the colour-colour analysis, we obtain a Mach number 
M = 4.58 + 1.09. 



The simulations performed by van Weeren et al. (201 lb i pro 



vide evidence that Mach numbers as high as 3 can also occur in 
merger shocks, but this is dependent on the model of the pre- 
shock temperature. The highly aligned magnetic fields mapped 
by van Weeren et al.| ( [2010[ ) also suggest compression of the ICM 
at the shock location. 

In terms of injection mechanisms, the models relate directly 
to physical processes operating at the shock front. The JP and KP 
models corresponds to a scenario at which electrons are accel- 
erated at the shock front and cool while streaming downstream. 
These two models have been described as having a 'single burst' 
of electrons which means that the particles in a mass element of 
the ICM are only accelerated while passing through the shock. 
The CI model implies that there is injection all the way down- 
stream. There are no published models for relic formation con- 
sistent with CI, but one could speculate that acceleration by 
small scale downstream turbulence induced by the shock (simi- 
lar to haloes) would have such a signature. The KGJP model can 
be reduced to a continuous sum of JP events with different age- 
ing. This is equivalent with the effect of mixing of different pop- 
ulations of electrons, which occurs because of projection effects. 
Because we see the relic side-on, there are parts of the spherical 
shock front that cause electron acceleration and are projected 
onto the downstream area. From our colour-colour plots, we ex- 
clude the continuous injection model as it does not match any of 
our data points. The data does not turn back to the line of equal 
low to high frequency spectral index, so it is also not well de- 
scribed by pure KP models. Despite the fact that the KP model 
does not fall within the error bars of most the points, we cannot 
fully rule it out. The SNR situation did not allow us to extend 
our spectral index coverage to steeper spectral index values in 



order to test whether the data turns back to the 



a]f 2 2 v The 



JP models do not provide a good fit of the extrema of the data. 
The KGJP model with a; n j = -0.6 provides the best description 
of the data. The deviations from this model can be explained by 
multiple factors. The effect of mixing is enhanced towards the 
centre of the cluster, where the steepest spectral indices appear. 
Another cause can be interfering emission from other sources, 
such as the radio source at the east end of the relic, the faint halo 
emission or unresolved compact sources. We also artificially in- 
troduced mixing when we binned pixels based on their spectral 
index. This acts as a complication to our analysis, since simple 
injection models are not applicable under these conditions. The 
additional mixing in the downstream area could cause the data 
to mimic a KP behaviour and turn back to the line of equal high 
and low frequency spectral index. 

2012[ ) produced simulated colour-colour plots 



Kang~et~aLl 



tailored to the physical prescriptions of C2242. While a direct 
comparison is difficult, since they used different reference fre- 
quencies, the general trends are the same. They also observe the 
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effects of radiative cooling curving the data away from the line 
of equal low to high frequency spectral index and notice that the 
trends in the colour-colour plots are highly dependent on projec- 
tion angle of the relic with respect to the sky. 

The results of the colour-colour analysis are similar to the 
ones obtained by van Weeren et al. (2012b} for cluster 1RXS 
J0603.3+421. This indicates that in both clusters the same un- 
derlying mechanisms play a role producing spectral shapes con- 
sistent with acceleration to a power-law and then ageing in the 
post-shock region (JP model) and mixing effects. 

4.2. Other relics 
4.2.1. Morphology 

By contrast to RN, the morphology of the counter relic RS 
changes significantly with frequency. At high radio frequencies, 
RS is connected with a fiat-spectrum tail (a = -1 .10±0.04, mea- 
sured between 1221 MHz and 1714 MHz) towards the south, 
but the tail is absent in the low frequency map. As mentioned in 
Sect. |3 . 1 . 1 1 the spectrum of the tail should be steeper than -1.4 
for it to be detected in the 153 and 323 MHz maps. While RS is 
separated from source J at 1714 MHz map, the two patches are 
connected into a single system at low frequencies. Figure A. 1 1 
reveals a point source with a red galaxy counterpart embedded 
in emission from J. The color of this galaxy is inconsistent with 
being a cluster member, but indicates it is a background source 
unrelated to the cluster. 

Source K towards the north of RS between sources G and 
F gets resolved into multiple compact, elliptical areas of emis- 
sion in the higher-resolution, low-frequency maps. There is also 
one head-tail source. All the components of source K have red 
optical counterparts. 

Rl is positioned at a similar projected physical distance from 
the cluster centre as source RN, suggesting that Rl and RN could 
have been accelerated by a common merger interaction. Rl also 
has an arc-like morphology, suggesting that the source is seen 
edge-on. By contrast, the patchy nature of R2 could result in 
a scenario where the emission is extended in the plane of the 
sky (i.e. seen face-on). Both sources do not have obvious optical 



counterparts coinciding with the radio peak (see Fig. A. 1 



In the 153 MHz map the two relics R2 and RN and source I 
become connected with diffuse emission which was interpreted 
as a radio halo ( |van Weeren et al.]|2010| l. R2 is most probably 
embedded in the halo that becomes prominent only at the lowest 
frequencies due to its steep spectrum. Diffuse source I could be 
a flatter-spectrum part of the halo. The measured fluxes within a 

wide box around the emission give a spectral index of 1 .25, 

a typical value for haloes. 

4.2.2. Spectral analysis 

The division between the southern relic and source J is clear ac- 
cording to their different spectral properties. Source J has a very 
steep spectral index, which explains why the source is barely no- 
ticeable in the high-frequency radio maps. RS and the southern 
radio tail seem to form a system of flatter emission. We speculate 
that RS and J is produced by two main broken shocks with dif- 
ferent acceleration parameters (i.e. Mach numbers). It is difficult 
to directly measure an injection spectral index for RS due to the 
disturbed nature of the relic. We can therefore, use the relation 
between the integrated (ai n j = — 1.29 ± 0.05) and injection spec- 
tral index to derive aj n j = -0.79 ± 0.05, which could be produced 
by a weaker shock than RN. 



Rl shows sign of spectral steepening from which supports 
the assumptions (see Sect. 4.2.1 1 that the source is viewed edge- 
on (-2.5 < a < -1.0). There is no clear spectral index gradient 
across R2. In a scenario where the relic is viewed close to face- 
on, this is naturally expected as we are probing though electron 
populations with different ages. The shock front may have also 
broken off into different fragments, as expected from simulations 
dPaul et al.|20TT) . 

4.2.3. Acceleration of emitting particles 

Apart from DSA/re-acceleration, there is one competing sce- 
nario for the formation of relics. Adiabatic compression requires 
the presence of a nearby radio source (which could be inac- 
tive now) to supply radio fossil plasma younger than 2 Gyr. 
Due to the old age of the plasma, the spectra of these so-called 
phoenixes are expected to be extremely curved and steep with 
a < -1.5 because of extensive energy losses. This scenario 
could explain why phoenixes do not trace an entire, spherical 
shock front, and predict most relics should be a few hundred kpc 



away from the cluster centres (Hoeft et al. 2004 Vazza et al. 
|20l2l ). 

RS, Rl, R2 could also trace shock waves but they are proba- 
bly produced by broken shock fronts and may also be viewed at 
an angle with respect to the line-of-sight which results in mixing 
of different electron populations and explains why clear spec- 
tral index trends are not visible. The prescriptions of adiabatic 
compression are not compatible with the observed position and 
integrated spectral index (~ —1.0) for these three relics. The 
flat integrated spectra of Rl and R2 (ai nt = -0.74,-0.90 re- 
spectively) are incompatible with pure (re-)acceleration theory, 
where the minimum injection index can be -0.5, given by an 
infinite Mach number. An explanation could be that Rl and R2 
are produced via oblique shocks, where the relationship between 



the injection and integrated spectral index does not hold (Kirk & 



Heavens 1989). As well, a flatter integrated spectrum could arise 



due spectral ageing occurring in inhomogeneous downstream ar- 
eas. A similar explanation has been given for the relic in cluster 
Abell 2256 ( |van Weeren et al.|20i2"c] l. It has to be noted that the 
scatter in the flux measurements is quite high, which makes the 
spectral index fits less reliable. One explanation could be that 
we are not measuring the total flux of the sources at all frequen- 
cies. Rl is the only relic that displays a possible spectral turnover 
which is particularly puzzling. For RS, we can use the relation- 
ship between the injection and integrated spectral index (derived 
in Sect. |4. 1 A\ , to derive a Mach number of M = 2.81+0.19. The 
diffuse patch of emission J has a very steep integrated spectrum 
(- 1 .53 ± 0.05) which suggests the source could be a phoenix em- 
bedded in the emission from RS. The high resolution 608 MHz 
image (Fig. [3} shows hints for a radio source with an optical 
counterpart (Fig. |A. 1 [ > that could supply the revived plasma. It 
is conceivable that source J and RS have been accelerated by a 
common shock front. 



4.3. Radio galaxies 
4.3.1. Morphology 

Sources B, C, D, E and F are clear head-tail radio galaxies which 
have bright elliptical optical counterparts, coinciding with the 
peak of the radio emission. A spectral index gradient is notice- 
able across all of them. The nuclei of the sources have a 0.5, 

whereas the tails steepen to values below -2.0. The most dra- 
matic steepening can be observed across source F, whose tail 
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reaches spectral index values of -3.0. Source G, which was be- 
low the noise level in the 2272 MHz map, has an ultra steep 
spectrum with a 1.8. This source has multiple possible opti- 
cal counterparts. The point source embedded in the western part 
of R2 is only properly resolved in the 323 and 608 MHz radio 
maps. Its spectral index is much flatter than the relic's and fairly 
constant across the source, which is consistent with point source 
properties. We could claim that source H is part of the northern 
relic, owing to its proximity. The source morphology in the high- 
resolution 608 MHz map provides evidence the source is a radio 
galaxy (either a disturbed tailed galaxy or a highly asymmet- 
ric, double-lobed source), which also has an optical counterpart 
visible in the Fig A. 1 In addition, its spectral properties are dif- 
ferent from the relic's, with no clear index gradient and flatter 
a ~ -0.75. 



4.3.2. Acceleration of emitting particles 

The tails of radio galaxies are indicative of their relative motion 
with respect to the ICM. Assuming a geometry in the plane of 
the sky ( |van Weeren et al.|2011b| l, we can directly interpret the 
bi-modality in the orientation of the tails as a tracer of the two 
merging clusters. The axes of sources B, C, D and E are aligned 
with the merger axis, with tails pointing away from the northern 
shock. Strong evidence for the merger scenario is also the fact 
that source F near RS, is directed the opposite way. In merging, 
highly disturbed galaxy clusters, after core passage, there is a 
displacement between the dark matter and the gas pertaining to 
each of the clusters. Dark matter and galaxies act like collision- 
less particles, while the IGM of the two merging clusters inter- 
acts through ram pressure and slows down. This means that the 
dark matter concentration can be found ahead of the ICM peak 
for each of the clusters. As well, the galaxies decouple from the 
plasma and run in front of the gas of their parent clusters. This 
effect has been observed in a number of objects, most famous of 
which being the Bullet cluster ( |Clowe e t al. 2006). Simulations 
by van Weeren et al. (2011b i and X-ray observations by Ogrean 
et al. 1 2.01 3) > indicate that C2242 is a post-core passage cluster 
merger, with the main cluster at the north and the subcluster at 
the south. Therefore, galaxies B, C, D and E are part of the main 
merging cluster, while galaxy F belongs to the infalling subclus- 
ter. The radio head-tail morphology together with the spectral 
index gradients provide additional support for the merger sce- 
nario. 



5. Conclusion 

We presented a combined low and high radio frequency GMRT 
and WSRT analysis of CIZA J2242.8+5301. All of our maps, 
present two counter relics, three smaller patches of diffuse emis- 
sion and several head-tail sources. The two main relics are 
perpendicular to the merging axis, while the radio sources are 
aligned with it. 

We show the first deep, low-frequency maps of this cluster 
and integrated spectra for the five relics and a number of radio 
galaxies. The spectrum of the northern relic is very well de- 
scribed by a power law with integrated spectral index of -1.06 ± 
0.05, which is consistent with the previous high frequency mea- 
surements. We derive a Mach number M = 4.58 ± 1.09. For the 
southern relic, we measure a steep spectral index of -1.29+0.05, 
which suggests the source has been accelerated by a weaker 
shock of M = 2.81 ± 0.19. Embedded in the emission from 
the southern relic, we discover a possible radio phoenix re- 
vived by adiabatic compression. This source has an extremely 



steep spectrum a mt = -1.53 + 0.04. The two smaller patches 
of diffuse emission have flat radio spectra (-0.74 + 0.07 and 
-0.90 ± 0.06), which suggests electrons are cooling in an in- 
homogeneous medium behind the shock. 

We find spectral index gradients across the northern and 
southern relic and the eastern patch of emission, indicative of 
shock acceleration with radiative losses into the downstream 
area. The smaller patch of radio emission to the south of the 
northern relic is probably seen face-on, as it does not have any 
clear spectral trends. 

All of the proposed models for the formation of relics predict 
spectral curvature in the downstream area. The faint nature of 
radio relics impose a low SNR working regime that has proven 
detection of spectral curvature extremely challenging. We are 
able to find, for the first time, a curvature gradient across a relic, 
which suggests active acceleration towards the northern edge of 
RN, while towards the cluster area, in the downstream region, 
the spectrum is dominated by radiation losses. This is confirmed 
by the colour-colour diagram, where we obtained similar results 
to |van Weeren et al.lpOTIb) for cluster 1RXS J0603.3+421. 

We also analysed the small scale structure along the relic and 
conclude that variations in the spectrum on scales of ~30 kpc are 
most probably caused by noise peaks in one of our radio maps. 
This warns against over-interpreting local structure in cluster 
relics. We speculate that the variations on scales on ~100 kpc 
along the length of the relic are to be trusted and are probably 
caused a more complex shock surface. 

We conclude that KGJP based models fits the data best, sug- 
gesting the electrons go through a phase of acceleration under 
a continuous isotropisation of the pitch angle to the magnetic 
field, followed by spectral ageing. Mixing of electron population 
is induced by projection and resolution effects. This also acts as 
a confirmation of the (re-)acceleration mechanism. 

Most of the radio galaxies in the cluster have a head-tail mor- 
phology and have clear spectral trends across their lobes. We 
discover a very interesting radio source with an extremely steep 

(a 1.8) and curved spectrum. We also uncover a bi-modality 

in the distribution of radio galaxies, supported by the trends in 
the spectral index and curvature in their tails. This provides ad- 
ditional support for the the merger scenario where the northern 
and southern relic are accelerated through the shock acceleration 
mechanism in two outward travelling shocks waves. 

The probing of the downstream area, where the radio spec- 
trum is curved and steepest, is limited by the noise levels cur- 
rently attainable. We therefore introduce a bias against faint 
emission by imposing an RMS cuf-off (e.g., Orru et aL||2007| l. 
In our particular case, the limiting factors were the intrinsic 
lower GMRT sensitivity at 153 MHz map and the 2272 MHz 
map, with poor SNR given the low surface brightness of the ex- 
tended emission at this frequency. In the future, telescopes like 
LOFAR and LWA, at the longest wavelength side, and the new 
upgraded Jansky Very Large Array, which samples the shortest 
radio wavelengths, are needed to better constrain the spectrum 
over more than three decades in frequencies. The high frequency 
data is crucial in determining which is the responsible injection 
mechanism, while the low frequency spectrum will further test 
shock acceleration up to the atmospheric cut-off frequency at 
10 MHz. A further test on the merging scenario could be a weak 
lensing or optical analysis of C2242 that would point out a pos- 
sible displacement of the mass concentration or galaxy distribu- 
tion, respectively, from the X-ray peak. 
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Appendix A: Optical counterparts of compact and 
diffuse sources 

To identify optical counterparts for the compact and diffuse ra- 
dio sources in the field, we overlaid 608 MHz contours on Isaac 
Newton Telescope (INT) optical images (Fig. 



A.l 



The images 

are composites of Wide-Field Camera images taken in the B, 
V, R and I filters between October 1-8, 2009. The total integra- 
tion times per filter is ~ 4000s. The data were flat-fielded and 
bias-corrected with IRAF ( Tody|19 93 ) and the mscred package 
( | Valde s 1998). The I and R band images were also fringe cor- 
rected. We removed cosmic rays and other artefacts by rejecting 
pixels above 3.0cr rms . 



Appendix B: Spectral index and curvature error 
maps 

In order to evaluate the uncertainties in the pixel-by-pixel spec- 
tral index and spectral curvature values, we produced error maps. 
Values in Figs. B.l and B.2 represent lcr uncertainties in the 



spectral index as resulting from the linear fit of the logarithm of 
the flux measurements as function of log frequency. 

In the fitting procedure, errors on the flux consisted of 10% 
absolute flux calibration added in quadrature to the RMS noise. 
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Fig.A.l. Optical images from the INT telescope. Contours from 608 MHz are overlaid at [4,8, 16,32,64, 128,256,512, 1024] x 
o"rms level. Sources are labelled as in Fig. |4] 

Figure [B3]presents the curvature uncertainty, which is the aver- 
age of the errors in the high and low frequency spectral index fits 
(these are subtracted to obtain the curvature values). 
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A. Stroe et al: Discovery of spectral curvature in the shock downstream region 
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Fig. B.l. Spectral index error map corresponding to Fig. 14 Fig. B.3. Spectral curvature error map corresponding to Fig. 16 
Figures are in the same order. 




Fig. B.2. Spectral index error map corresponding to Fig. 15 
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